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Abstract 

We present the numerical solution of the full gap equation in a weak coupling 
constant g. It is found that the standard approximations to derive the gap 
equation to the leading order of coupling constant are essential for a secure 
numerical evaluation of the logarithmic singularity with a small coupling con- 
stant. The approximate integral gap equation with a very small g should be 
inverted to a soft integral equation to smooth the logarithmic singularity near 
the Fermi surface. The full gap equation is solved for a rather large coupling 
constant g > 2.0. The approximate and soft integral gap equations are solved 
for small g values. When their solutions are extrapolated to larger g values, 
they coincide the full gap equation solution near the Fermi surface. Further- 
more, the analytical solution matches the numerical one up to the order one 
0(1). Our results confirm the previous estimates that the gap energy is of 
the order tens to 100 MeV for the chemical potential fi < 1000 MeV. They 
also support the validity of leading approximations applied to the full gap 
equation to derive the soft integral gap equation and its analytical solution 
near the Fermi surface. 
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I. INTRODUCTION 



The deconfined quark matter is expected at baryon densities much larger than the normal 
nuclear matter density. Any attractive quark-quark interaction causes instability of the 
Fermi surface by formation of Cooper pairs and leads to the color superconductivity phase 
|l|J^Q] in a very dense quark matter at low temperatures. However, the quark matter at 
relatively low densities may have properties qualitatively different from those for BCS-type 
superconductor. It is argued that the nuclear matter at low density might be continuously 
connected to the quark matter at high density without any phase transition. The calculation 
of superconducting gap is important to study many properties of hadronic matter at high 
densities 

In QCD, the single-gluon exchange between quarks of different color generates an at- 
tractive interaction in the color antitriplet channel. The scattering through a single-gluon 
exchange strongly correlates the direction of the in- and outgoing quarks. However, there is 
a logarithmic divergence for the forward angle scattering. This implies that the gap equation 
is not an exponential in 1/(7^ Q, as in BCS like theories, but only in l/^f where g is the 
coupling constant. It is found for the nonzero quark density at low temperature, the free 
energy is an expansion in g"^ g), and appears well behaved even for much larger coupling 
constant < 4 |]^. 

The earlier estimations of the color superconductivity gap were found rather small A ^ 
(few) MeV for quark chemical potentials that are of physical interest /i < lOOOMel^ 
Recently, it is shown that the gap is large and of order 100 MeV for yU < 1000 MeV 
These large estimations may lead to observable effects in the core of neutron stars and 
the physics of heavy ion collisions P,|5|,p!0[| . The reliable estimation of gap equation plays 
an essential role in the phenomenological calculations involve the color superconductivity 
physics. 

The gap equation has been studied numerically for a quantitatively estimation for a 
wide range of the chemical potential in the time-like energy space as well as momentum 



space piip!2| , p!3| . In these studies, the gluon propagator has been specified by Thomas- Fermi 
screening and Landau damping defined by the hard dense loop approximation (HDL). Fur- 
thermore, the resultant equation has been simplified by approximating the angular inte- 
gration near the Fermi surface |T^. Nonetheless, the gap equation is highly singular near 
the Fermi surface due to the correlation of in- and outgoing quarks direction. Therefore, 
any approximation may affect significantly the gap value in particular for a moderate quark 
chemical potential. In the present work, we shall adopt the spectral representation of the 
gluon propagator |I^Jl5[] and solve the gap equation exactly numerically without any further 
approximation. The spectral representation is more rigorous than previous approximation 
based on screening gluon propagator. 

Recently, Pisarski and Rischke [|I^,|T^,|T^,|I^] have derived the full gap equation in the 
momentum space for two massless quark flavors with a total spin J = at an arbitrary 
temperature T. They adopted the gluon propagator in the spectral representation derived 
from HDL approximation. They derived the approximate gap equation to the leading order 
of coupling constant g. It is given by Eq.(72) in Ref. fl^. This equation has a logarithm 



dependence on the quark energy under the integral. Furthermore, they have approximated 
the logarithmic term under the integral to derive a soft integral equation that is given by 
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Eq.(80) in Ref. [0. Hereafter, the approximated integral gap equation is called Eq. (I) 
while the soft integral gap equation is called Eq. (II). Finally, they have inverted the soft 
integral equation to an ordinary differential equation to derive the analytical solution. 

In the present work, the full gap equation is solved numerically for a rather large coupling 
constant g = 2 — A which corresponds to the chemical potential fi < 1000 MeV. We confirm 
the previous estimates that for the chemical potential of order 1 GeV, the gap amplitude 
can be of order tens to 100 MeV. We also study the numerical solution for a small coupling 
constant g < 2.5 by solving Eq. (I) as well as Eq. (II). We compare the numerical solutions 
with the analytical one to the leading order logarithmic accuracy. This paper is organized 
as follows. In Sec. II, we present the full gap equation and its approximated versions. The 
results and the conclusions are presented in Sec. III. 



II. THE GAP EQUATION 



The full gap equation for the color superconductivity in the weak coupling with total 
spin zero and two massless flavors in a dense QCD reads ||T^ 



m = -/ 
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where 



{q — + 10(g) p. The spectral densities, pi and pt, are given by 



(6) 



The explicit expressions for the spectral densities pP°^^'^^^ qj-q given by Eqns(34b-e) in Ref. 



T^ . The co'i,t(p) are the solutions of the dispersion relations for longitudinal and transverse 



gluons which are given by Eqns(37a-b) in Ref. |]T^ . 

In HDL approximation there are three scales in cold dense QCD, the chemical potential 



/i, the gluon mass nig ~ gp and the color superconductivity condensate 
and they are naturally ordered p» nig » 



~ /iexp(-l/^) 
Pisarski and Rischke [ITli have derived the 



approximate integral gap equation by using the assumptions of these scales from the full 
gap equation. This approximate integral gap equation is derived in the momentum space 
unlike to the Eliashberg equation given by Schafer and Wilczak is derived in the time-like 
energy space. 

The approximate integral equation, Eq.(I), reads 



m 
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This equation can be simplified to a soft integral equation by 



replacing the logarithm under the integral by 
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In 1^^^ e{q-k)+\Yi (^^^ e{k-q). 



Hence, the soft integral gap equation, Eq.(II), becomes 
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(9) 



where g 

The soft integral gap equation, Eq.(H), is inverted to an ordinary differential equation 
which has a simple solution to the leading order in 'g and Fermi surface x* = vr/2^. The 
solution for zero temperature case reads 



26/iexp(-7r/(2^))sin(^a;) x C(l) 



(10) 



where x 



ln( 



2bii 



Furthermore, the analytical solution of the gap equation in the 



time-like energy space, Eliashberg equation, is given by 0] 



IT 



bpexp{-—) X 0(1) 



(11) 



where the overall coefficient is correct up to a prefactor of order one "(9(1)". Since the ana- 
lytical solution of the gap equation to leading order of coupling constant g is undetermined 
by a prefactor of order one, we assume that the analytical solution at the Fermi surface 
reads 



TT 
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(12) 



where ( is constant of order one. 
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III. RESULTS AND CONCLUSIONS 



We have studied the numerical solution of the full gap equation and the approximations 
validity for Eqns. (1) and (11) near the Fermi surface. The full gap equation is solved in 
momentum space for a rather large coupling constant g > 2.0. The approximate versions 
of the full gap equation (i.e. Eq.(I) and (2)) arc solved for small coupling constant values 
g < 1. Furthermore, we have assumed that Eqns. (1) and (11) are extrapolated to large 
coupling constant values and then we have presented their numerical solutions for coupling 
constant values 1.0 < g < 3.0. 

Fig. 1 displays the relative gap amplitude 0(|k|)/0o versus the scaled Fermi momentum 
k/iJ,. It is shown that the solution of Eq. (11) have a sharp peak near the Fermi surface 
for small coupling constant values g < 1. The peak around the Fermi surface for Eq.(ll) 
is much more sharp than that for Eq.(I) with a small coupling constant g < 1. Hence, the 
replacement of logarithm under the integral in Eq. (1) by step functions in Eq. (11) secures 
the singularity problem more thoroughly than the Cauchy numerical integration in Eq. (I) for 
a small coupling constant g < 1.0. For the large coupling constant, g > 1, both equations 
(1) and (11) give similar results. As we shall show below, Eq. (1) gives much a smaller 
amplitude than that for Eq. (II) with a small coupling constant g < I. The reason for 
that is simple. The gap equation appears very singular near the Fermi surface for a small 
coupling constant. Hence the approximation in Eq. (II) evaluates the gap amplitude value 
near the Fermi surface more thoroughly than numerical Cauchy integration. Furthermore, 
the gap amplitude in Eq. (11) approaches the same behavior for that in Eq. (1) with large 
coupling constants g > I. Therefore, the numerical Cauchy integration in Eq. (I) becomes 
unreliable for a very small coupling constant g < 1. 

Fig. 2 displays the gap amphtude 0/// for the full gap equation versus the relative 
momentum k/n with rather large coupling constant values g > 2.0 that correspond chemical 
potentials of the order 1 GeV. It is shown that the gap amplitude has a peak near Fermi 
surface and the function becomes more smooth as the coupling constant g increases. 

Fig. 3 displays the gap amplitude versus the logarithmic parameter x for the gap Eqs. 
(I) and (II) as well as the analytical solution normahzed to the amphtude of Eq. (II). 
The maximum value represents the gap amplitude near Fermi surface. Hence the values of 
4>{xo) I /i and xq arc important near Fermi surface k = fi. It is shown that for small coupling 
constant values g < 1, the analytical solution and numerical solution of Eq. (II) have almost 
the same behavior while the numerical solution of Eq. (I) has a significant deviation. As 
the coupling constant g increases, the numerical solution of (I) and (II) approaches each 
other and almost coincides for g > 1.5. However, the analytical solution deviates from the 
numerical one (II) significantly as the coupling constant exceeds g > 1.5. This means that 
Eq. (I) and Eq. (II) becomes almost the same for a large g and the standard approximation 
which are applied to Eq. (II) to derive the analytical solution fails at large coupling constant. 

In Fig. 4, the gap amplitude 0(x)//x versus the logarithmic parameter x with coupling 
constant g = 2.5 is displayed for the full gap equation as well as for Eqns. (I) and (II). 
It is shown that the amplitude of the full gap equation approaches the values of Eqns. (I) 
and (2) more closely near the Fermi surface (i.e. at the hill of the curves) while it deviates 
significantly for momentum far away from the Fermi surface. Eqns. (I) and (II) solutions 
have almost similar behavior for the most values oi x. It is interesting to note here that the 
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validity of Eqns. (I) and (II) is restricted only near the Fermi surface and any solution far 
away from the Fermi surface isn't guaranteed. Hence, Eqns. (I) and (II) give asymptotically 
accurate solutions near the Fermi surface for large coupling constant values g >2. However, 
Eq. (I) solution seems accurate for a moderate couphng constant g since the singularity 
near the Fermi surface becomes less severe. In the other hand, the solution of Eq. (II) seems 
more accurate than that for Eq. (I) for small coupling constant since it seems to account 
the singularity near Fermi surface more precisely. 

Fig. 5 displays the gap amplitude (po/ a.t the Fermi surface versus the coupling constant 
g for the full gap equation and its approximation to Eqns. (I) and (II) as well as the analytical 
solution with prefactors ( = 0.5 and 1.0. We have shown the solutions for Eqns. (I) and 
(II) with g > 0.7 and the solution of the full gap equation with g > 2.0. The sharp peak 
near the Fermi surface for Eq. (II) is larger than that for Eq. (I) for a coupling constant 
g < 1 since Eq. (II) is expected to secure the singularity near the Fermi surface better than 
Eq. (I) as mentioned above. However, the solutions (I) and (II) become close to each other 
for a wide range of coupling constant values 1 < g < 3. The solutions of Eqns. (I) and 
(II) match the solution of full gap equation for 2 < g < 3.0. They deviate from the full 
gap equation for g > 3.0. The full gap equation solution, for g > 3.0, starts to increase 
linear logarithmic with respect to g. We also display the analytical solution with prefactor 
coefficients ( = 0.5 and 2.0. The analytical solution with ( = 2.0 seems to fit the solution 
of Eq. (I) for small g < 1. When the coupling constant g increases, the analytical solution 
with a prefactor coefficient zeta = 1/2 seems to fit the solutions of Eqns. (I) and (II) as well 
as the full gap equation for a range 1 < g < 3.0. However, the analytical solution fails to 
fit the full gap one with a specific prefactor coefficient for the large coupling constant range 
g > 3.0. Therefore, the analytical solution approximation breaks down for a rather large 
coupling constant g > 3.0. 

Fig. 6 displays the amplitude of full gap equation at the Fermi surface versus the 
chemical potential /i in the range 550 MeV < fi < 800 MeV. We have used the running 
coupling constant defined by g"^ = (87r^//5o) x 1/ ln(/i^/AQC'D) where jSo = {llNc — 2Nf)/3. 
The QCD scale is taken as Aqcd = 400 MeV. The value of the gap amplitude is found of 
the order of tens to 100 MeV with quark chemical potential that may exist in neutron stars. 
This large value of the gap amplitude near the Fermi surface confirms the earlier estimates 

0- 

In summary, we have studied the numerically solution of the full gap equation and its 
approximations. The numerical results are compared to the analytical one with the leading 
order g approximation. It is found that the gap equation has a sharp peak near the Fermi 
surface for a small coupling constant g < 1. The singularity at the Fermi surface is hard to 
secure numerically in the full gap equation with small g. The approximation to a leading 
coupling constant g near the Fermi surface is essential to smear and evaluate this singularity 
in Eqns. (I) and (II). When Eqns. (I) and (II) are extrapolated to large coupling constant 
values g > 2, they match the results of the full gap equation near the Fermi surface. It is 
shown that the analytical solution coincides the solution of Eq. (II) for small g. However, 
this analytical solution approaches those of Eq. (I) and full gap equation with a rather 
large coupling constant g < 1.5 only near Fermi surface. This justifies the validity of 
approximation hold in Eqns. (I) and (II) near Fermi surface. 
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FIG. 1. The color superconductivity gap amplitude (})/(pQ versus k/^ with various coupling 
constant g. There is a peak at the Fermi surface momentum k/jj, = 1. (a) Eq. (I), (b) Eq. (II). 
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FIG. 2. The color superconductivity gap amplitude ^o/a* versus k/ (i for the full gap equation 
with various coupling constant g. 
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FIG. 3. The gap amplitude <p{x) as function of the variable x. The solid line is a solution of 
Eq. (I) and the dashed line is a solution of Eq. (II) . The dot-dashed line is an analytic solution to 
the leading order of coupling constant g. 
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FIG. 4. The gap amplitude <^(x) as function of the variable x with coupling constant g = 2.5. 
The solid line is a solution of the full gap equation. The dashed line is a solution of Eq. (I) as well 
as the doted line is a solution of Eq. (II) . 
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FIG. 5. The color superconductivity gap amplitude </>o//U versus the coupling constant g using 
several approximations. The solid line is a solution of Eq. (I) and the dashed line is a solution 
of Eq. (II) while the dotted-dashed line is a solution of the full gap equation. The light solid 
and dotted lines are the analytic solutions to leading order of coupling constant g multiplied by 
prefactor coefficient ( = 0.5 and ^ = 2.0, respectively. 
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FIG. 6. The color superconductivity gap amplitude versus the chemical potential for 
the full gap equation. The coupling constant g is identified with a standard logarithmic running 
couphng with QCD scale I^qcd = 400 MeV. 
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